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ABSTRACT 

We present Hubble ACS images of thirteen dust reddened Type-1 quasars selected from the 
FIRST/2MASS Red Quasar Survey. These quasars have high intrinsic luminosities after correction 
for dust obscuration (—23.5 > Mb > —26.2 from K-magnitude) . The images show strong evidence 
of recent or ongoing interaction in eleven of the thirteen cases, even before the quasar nucleus is 
subtracted. None of the host galaxies are well fit by a simple elliptical profile. The fraction of quasars 
showing interaction is significantly higher than the 30% seen in samples of host galaxies of normal, 
unobscured quasars. There is a weak correlation between the amount of dust reddening and the mag- 
nitude of interaction in the host galaxy, measured using the Gini coefficient and the Concentration 
index. Although few host galaxy studies of normal quasars are matched to ours in intrinsic quasar 
luminosity, no evidence has been found for a strong dependence of merger activity on host luminosity 
in samples of the host galaxies of normal quasars. We thus believe that the high merger fraction in 
our sample is related to their obscured nature, with a significant amount of reddening occurring in 
the host galaxy. The red quasar phenomenon seems to have an evolutionary explanation, with the 
young quasar spending the early part of its lifetime enshrouded in an interacting galaxy. This might 
be further indication of a link between AGN and starburst galaxies. 
Subject headings: quasars: general — galaxies: active, evolution, interactions 



1. INTRODUCTION 

The existence of a link between AGN and star- 
bursts/Ultraluminous Infrared Galaxies (ULIRGs) in 
the local Universe has b e en d is cussed for a long 
time (e.g. 'Heckman ct al.' (1989"); ISanders fc Mirabel 
(|1996D : [cTd Fcrnandcs ct al. (2001)). Further evidence 
for such a link comes from the tight correlation be- 
tween nuclear black hole mass and b ulge mass (e.g 
iMagorrian et al.l (|1998D : iTremaine et all 120021) and the 
black hole mass - host galaxy velocity dispersi on (M- 
a) rel a tion in local ga l axies (e.g. [Fcrrarcse & MerrittI 
(|2000( l: lGebhardt et al.l (|200G[ )). which suggest that star- 
formation and accretion onto black holes in galax- 
ies are intimately connected. However, HST imaging 
of the host galaxies of luminous quasars to moderate 
depths showed that most of them have the morphol- 
ogy of elliptical gala .xies, with little evi dence for ongo- 
ing star formation (jDunlop et al.l 120031 ) . At low red- 
shift only 30% of the host galaxie s show obvious signs 
of disturbance (jCuvon et al.l 120061 ) . Although quasars 
with infra red excesses are thought to host starbursts 
(jCanalizo fc Stocktonll2001h . these are a small minority 
of the optically-selected quasar population. 

Recent surveys have shown that optically-selected 
quasars comprise less t han half of the total quasar popu - 
lation (Martfncz-Sansi gre et al.ll2005Hs"tern et al.ll2005f ). 
Optical quasar surveys tend to miss dust-reddened 
quasa rs which have been found either wi th infrared sur- 
veys (jCutri et al.ll2001t iLacv et al.ll2004f l. radio surveys 
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(I White et al.l[200l . hard X -ray surveys or man et al.l 
12002) or surveys for high-ionization narrow-line objects 
(Zakamska ct al. 2006). Therefore they represent a new 
and largely uninvestigated quasar population, which may 
have many members that are at an earlier stage in their 
quasar activity in which dust and gas debris from the 
merger block the view of the central AGN. This fits with 
population synthesis models of the X-ray background 
which predict that a large fraction (close to 80%) of 
the accretion process in AGN is obscured and only the 
most energetic photoi is (> 20 keV) can penetrate the ob- 
scuring dust barrier fCiUi ct al.' ' 20011 : llJeda et al.ll2003t 
GiUi. Comastri & Hasinger 2007)). 

Dusty quasars are generally divided into Type-2 
quasars, which show only narrow- line emission in the 
rest-frame optical, and have inferred extinctions towards 
the nucleus of Ay ~ 5 — 100, and lightly dust-reddened 
Type-1 quasars, which still show broad lines and quasar 
continua in the optical, but have significant reddening 
relative to a normal quasars {Ay ~ 1 — 5). 

Surveys in the mid-infrared have been able to make 
approximate c ensus of th e obscured quasar population. 
iBrow n et all (|2006l ) and iLacv et al.l (|2007l ) find that 
the fraction of lightly dust-reddened Type-1 quasars is 
between 20 and 30% of the total Type-1 population. 
iLacv et aD (|2007| ) also shows that the Type-1 population 
(unreddened plus lightly reddened) is about 50% of the 
total quasar population, the remainder being Type-2s. 

The nature of both populations remains uncertain. At 
present, there are two competing models to explain the 
obscuration in dusty quasars: orientation vs. evolution. 
In the orientation model, dust reddened quasars are ob- 
serv ed side-on to a d isk or torus of obscuring material 
(e.g lAntonuccil (|1993i )). whereas in the evolution model, 
dusty quasars are viewed as an early stage in the life 
of a quasar, where the nucleus is obscured by debris 
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Fig. 1. — Spectra of the FIRST/2MASS red quasars observed with HST. The black Hne in the quasar spectrum, the blue line the FBQS 
quasar composite redshifted to the quasar's redshift and the red line after applying a SMC reddening law to the FBQS composite. 



from a merger and starburst which triggers the nuclear 
activity, prior to qu asar winds expeUing the dust (e.g. 
ISanders et al.l ()1989f )). Of particular interest are the 
class of lightly reddened quasars, where the reddening 
is smaller than that expected from a torus seen edge-on. 
These objects are good candidates for reddening by dust 
in the host galaxy, but may also correspond to objects in 
which the line of sight grazes the edge of the torus. 

Recent simulations of galaxy mergers tend to sup- 
port the evolution model, and suggest that the active 
nucleus is obscured for a long time before the feed- 
back from the accret ion disperses the obscuring material 
([Hopkins et al.l[2005h . In this picture, the quasar is seen 
in the optical regime only at the end of its lifetime, when 
it is the most luminous, radiating close to the Eddington 
limit and therefore being capable of expelling the dust 
([Hopkins et al.l[2006h . A large population of these young, 
obscured, underluminous quasars would also account for 
the mostly hard X-ray background; however there is still 
debate on the normalization factor of the hard X-ray 



background, that is the fraction of these sources is still 
unknown. Observations of Type-2 quasars discovered in 
the Sloan Digital Sky Survey (SDSS), however, tend to 
support the orientation hypothesis, with broad lines seen 
in polarized light, and host galax ies consisten t with those 
of normal Type-1 quasars ( Zakamska et al.l 12006.) . Sim- 
ilarly, an HST study of lightly-reddened Type-1 quasars 
selected from 2MASS found no significant difference in 
the properties of host galaxies of infrared sel ected quasars 
from quasars selected from other methods fM arble et al.l 
l2003f ). though ground-based studi es of a similar sample- 
found more signs of interactions ([Hutchings et al.ll2003l 

This paper focuses on HST observations of a sample of 
lightly dust-reddened Type-1 quasars with 0.4 < z < 1, 
selected using a combination of the FIRST radio survey 
and the 2MASS infrared survey. The primary objective 
is to study the host galaxies of red quasars and to assess 
if their morphologies show signs of recent interactions. 
We adopt a flat universe, Hq = 70 km s~^ Mpc~^, Ha = 
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TABLE 1 

Source properties 



Radio Flux E(_B — V) Luminosity 



Source 


z 


K-mag 


1.4GHz (mjy) 


(rest' 




Spectrum 


/ K-mag 


F2M0729+3336 


0.954 


14.5 


3.3 


0.83 ± 


0.22 


-25.0 


/ 


-25.6 


F2M0825+4716 


0.804 


14.1 


61.1 


0.50 ± 


0.35 


-21.9 


/ 


-25.8 


F2M0830+3759 


0.413 


14.6 


6.4 


0.80 ± 


0.15 


-22.7 


/ 


-23.5 


F2M0834+3506 


0.470 


14.7 


1.2 


0.53 ± 


0.10 


-22.9 


/ 


-23.7 


F2M0841+3604 


0.552 


14.9 


6.5 


1.05 ± 


0.64 


-21.0 


/ 


-23.9 


F2M0915+2418 


0.842 


13.8 


9.8 


0.53 ± 


0.36 


-23.7 


/ 


-26.2 


F2M1012+2825 


0.937 


15.2 


9.2 


0.66 ± 


0.12 


-24.4 


/ 


-25.0 


F2M1113+1244 


0.681 


13.7 


3.0 


1.01 ± 


0.24 


-25.1 


/ 


-25.5 


F2M1118-0033 


0.686 


14.6 


1.3 


0.72 ± 


0.26 


-23.0 


/ 


-24.6 


F2M1151+5359 


0.780 


15.1 


3.5 


0.47 ± 


0.13 


-23.3 


/ 


-24.7 


F2M1507+3129 


0.988 


15.1 


7.8 


0.54 ± 


0.13 


-24.3 


/ 


-24.8 


F2M1532+2415 


0.564 


15.0 


7.4 


0.90 ± 


0.53 


-20.5 


/ 


-23.9 


F2M1656+3821 


0.732 


15.1 


4.1 


0.55 ± 


0.57 


-20.9 


/ 


-24.5 



0.7 cosmology. 

2. BACKGROUND AND SAMPLE SELECTION 

2.1. The F2M sample 

We have selected a sample of luminous, dust-reddened 
quas ars using a com bination of 2MASS infrared sur- 
vey (iSkrutskie et a.1.1 [l995.) and the FIRST radio sur- 
vey (|Becker et alj |1995[) in order to answer some of 
the questions posed above. The bulk of these dust- 
reddened Type-1 quasars have reddenings around E{B — 
^) ^ 1-5, which makes them much more heavily red- 
dened than other red quasar samples (jCutri et al.ll2001l : 
[Richards et al.|[2003l ). but they do not show the amount 
of obscuration as true Type-2 quasars. The method- 
ology and preliminary re s ults o f this sample are de- 
sc ribed in 'Glikman et al.l (|2004| ) and the final results 
in I Glikman et al. (2003). We have used two different 
color selection techniques, both of which seem to be 
highly effective at finding red quasars {J — K > 1.7 and 
i? — -ft' > 4; or i? — > 5). We have followed up our 
candidates with Keck (ESI) and IRTF (Spex) to obtain 
spectra in the optical and near-infrared. Some of these 
objects show only narrow lines in the observed-frame op- 
tical, but spectroscopy of the quasars in the near-IR have 
show very broad Ha or Paschen lines. So far, about 100 
red quasars have been found with this method, which 
woul d have been missed by traditional optical quasar sur- 
veys (jGlikman et al.|[2007t ). Following their conventions, 
objects from the FIRST-2MASS red quasars survey are 
named F2M. 

In contras t to the 2MASS quasar sample of 
iMarble et ed] jlOOl), which were mostly low redshift, 
low luminosity and selected by their near-infrared col- 
ors alone, the red quasars in our FIRST/2MASS study 
have a median redshift of 0.7. We therefore expect the 
F2M sample to have have higher luminosity and perhaps 
more star-formation in the host galaxies than other red 
quasar samples at lower redshift. 

While the selection of red quasars based on R — K 
colors could include objects which are red due to galaxy 
starlight in the infrared, the objects in this sample are 
intrinsic high luminosity quasars, in which the galaxy 
light should add a negligible contribution in the K-band. 
This becomes evident when we measure the reddening of 
the quasars. 



2.2. Sample spectra 

From the F2M sample, we chose a subsample of 13 
0.4 < z < 1.0 red quasars to image to inspect their host 
galaxy properties. In this section we will describe their 
optical spectra, taken mostly with ESI on Keck. 

We fit the optical red quasar spectra to a model of an 
unreddened q uasar. For that we use the FBQS compos- 
ite spectrum (iBrotherton et al.l 2001 ) with a SMC dust 
reddening law from lFitzpatrickl ( 199S ) using the relation. 



Fo(A) =F,(A)10^(^- 



y)fe(A) 



(1) 



with Fc{X) the composite (FBQS) and Fo(A) the 
observed spectrum. Typically the SMC extinction 
curve fc(A) lacks a significant 2175 A bump seen 
in Galactic dust. We fo llow the conventions from 
iFitzpatrick fc Massal ()1990( ) to obtain the extinction 
curve k{X), with Ry — 3.1 (Galactic value). 

Figure [1] shows the spectra of our quasars. The black 
line is the observed spectrum, the blue line is the FBQS 
composite redshifted and normalized to the quasar flux 
and the red line is the FBQS composite with the deduced 
reddening correction applied. The derived E{B — V) are 
shown in Table [U the values have such high errors in 
E{B — V) either because there was low signal to noise 
in the spectrum or the rms of the fit was not very good. 
Because there can be significant host galaxy contribution 
to the optical spectra, especially in the UV, our E(B — 
V) est imates can serve as a lower limit. iGlikman et al.l 
(|2007f ) also includes the near-infrared SPEX spectra to 
their fitting, so in some cases there is some discrepancy 
between their quoted values and ours. Their deduction 
of the reddening is beyond the scope of this paper, so we 
will only take the optical spectrum for an estimate of the 
reddening of the total system. 

For most of the quasars the reddening fits are con- 
sistent with dust-reddening and very little red starlight 
from the host galaxy is required. However, this is 
not the case for F2M0825-I-4716, F2M0915-I-2418 and 
F2M1151-h5359 in which the the fit for E{B - V) broke 
down. These objects will be important later when we 
inspect the properties of their host galaxies. Note also 
that the reddenings are significant enough that the red 
colors of the quasars could not solely have come from a 
"red" spectral slope. Even though the F2M sample is 
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Fig. 2. — Color composites of the 13 red quasars. We use the I-band data for the red part and the g-band for the blue and green parts of 
the composite. The images are 7 X 7 . The g -band cuts were set so that it the blue/green part of the image is somewhat more enhanced 
in comparison to the 7c-band data; the actual details from the host galaxies can be seen. The true color images would appear redder than 
shown here. For a quasar at redshift of z ~ 0.7, which is typical for our sample, 1 arcsecond represents about 7 kpc. 



radio selected, the radio fluxes of the objects (Table [T]) 
are also weak; the quasars are radio intermediate at best, 
so a strong synchrotron component showing up in the in- 
frared is unlikely to have much contribution to the red 
color. 

Once we have the reddening E{B — V), we can also 
derive the intrinsic luminosities of the F2M quasars from 
the spectrum (corrected for obscuration). We also cal- 
culated the luminosities of the quasars from the K-band 
magnitude assuming that in the K-band there is no host 



galaxy contribution. The range between the absolute 
magnitude derived from the spectrum and from the K- 
band magnitude sometimes is quite large implying that 
there either were some large slit losses or that even in the 
K-band there is significant host galaxy contribution, see 
Table[TJ We are more inclined to use the derived luminos- 
ity from the K-band, since the slit will miss the obscured 
AGN often and as mentioned before, the E{B~V) values 
are only a lower limit as there is host galaxy contribu- 
tion to the spectrum and we don't have the full wave- 
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Fig. 2. — cont. 



length coverage ranging into the Mid-Infrared to get the 
true continuum shape. The fact that the K-Band mag- 
nitude is so bright imphcs that the reddenings must be 
larger. What we can say is that the quasars are more 
luminous than usual, having a luminosity range of —23.5 
> Mb > —26.2 derived from the K-Band magnitude. 
Typical quasar luminosity functions at similar redshifts 
as our sample have only the highest luminosity quasars 
above Mb ~ —24.5 (our mean value), with a steepening 
of the function at arou nd Mr ^ ~^3-5 for q uasars at a 
redshift of z - 0.7 (e.g. iHopkins et all (|2007t )). The high 
values of Mb we find for our quasars suggests that we 
are missing the majority of this type of objects. Thus, 
because of the shallowness of the 2MASS survey we are 
only able to probe the tip of the "red quasar iceberg" . 

3. HUBBLE ACS OBSERVATIONS 

We followed up this subsample of 13 F2M red quasars 
with the ACS Wide Field Camera on HST (GO-10412). 
They were imaged with the F475W and the F814W fil- 
ters, which roughly correspond to g and Ic filters. The 
redshift range was chosen such that all our objects had 
their intrinsic luminosities well above the quasar/Seyfert 
divide (so that the near-IR magnitudes used in their se- 
lection would not have much contribution from the host 
galaxies) , and all were at low enough redshift that emis- 



sion above rest-frame 4000 A fell within the bandpass 
of the F814W filter on ACS. The observing dates and 
corresponding exposure times are presented in Table [21 

Most of the objects were observed for one HST or- 
bit. However, as the sample of red quasars spanned a 
significant range of redshift, two high redshift quasars 
(z > 0.9) were observed for two orbits each to ensure a 
more uniform surface brightness limit for the sample as a 
whole. These images were drizzled with the python pro- 
gram multidrizzle in the usual manner. For the images 
with only one orbit we used the readily available driz- 
zled products from the archive, since they did not have 
significant cosmic ray contamination. 

Color composite images of the red quasars are shown in 
Figure[2l Inspecting the images by eye, one can see much 
more information of the host galaxy than one would from 
a blue quasar in which the quasar light dominates. How- 
ever, the images reveal much more interaction than the 
usual fraction of 30% (jMarble et al.ll2003t iGuvon et all 
l2006( l with tidal tails and many irregularities present in 
the host galaxies. 

Some of the images show several compact knots or nu- 
clei. If the merger was driven by two galaxies with central 
black holes, it is possible that the two of them ignite as 
AGN, before merging into a massive active galaxy with 
only one black hole in the center. This has been ob- 




TABLE 2 
Journal of HST observations 











Exposure 


Time 


Source 


Obs. 


Date 


red (s) / 


blue (s) 


F2M0729+3336 


Oct 


24 


2005 


4650 


/ 


4760 


F2M0825+4716 


Oct 


19 


2005 


1820 


/ 


1924 


F2M0830+3759 


Apr 


18 


2005 


1720 


/ 


1832 


F2M0834+3506 


Apr 


18 


2005 


1620 


/ 


1720 


F2M0841+3604 


Apr 


24 


2005 


1600 


/ 


1716 


F2M0915+2418 


Jun 


05 


2005 


1672 


/ 


1776 


F2M1012+2825 


Jun 


13 


2005 


2190 


/ 


2250 


F2M1113+1244 


May 


26 


2005 


1660 


/ 


1760 


F2M1118-0033 


Jul 


14 


2005 


1640 


/ 


1752 


F2M1151+5359 


Apr 


20 


2005 


1860 


/ 


1972 


F2M1507+3129 


Jul 


28 


2005 


2290 


/ 


2340 


F2M1532+2415 


Apr 


24 


2005 


1736 


/ 


1840 


F2M1656+3821 


Apr 


19 


2005 


1760 


/ 


1876 



served in NGC6240, an Ultraluminous I nfrared Galaxy 
at very low redshift ()Komossa et al.|[2003D . Table |4] sum- 
marizes the properties of the quasars that show more 
than one candidate nucleus. F2M0841-I-3604, which has 
the largest angular separation, was observed with Chan- 
dra, but only had 7 (very hard) counts, which is not 
enough to resolve both components. Two of the objects 
also have double-peaked broad lines, F2M0825-I-4716 and 
F2M1507-I-3129. Such emission lines are not uncommon 



in broad-line radio galaxies, and ar e thought to be due to 
a disk -like BLR seen edge-on (e.g. lEracleous &: HalpernI 
(|1994f )). however, they could also be from two separate 
active nuclei close to merger. Unfortunately, neither of 
these show multiple nuclei in their images. 

We then performed photometry measurements on the 
ima ges in the two bands. For this we used Sextrac- 
tor (iBertin fc ArnoutillQQGT) w ith parameters set as de- 
scribed in lBenitez et al.l ( 20041 ). We used the AB mag- 
nitude zeropoints from the ACS webpage. We extracted 
two magnitudes this way: first the isophotal corrected 
magnitude (assuming a symmetric Gaussian profile for 
the object) and then a 3" aperture magnitude. A large 
difference between the two may indicate that a lot of the 
light comes from a galaxy component outside the nu- 
cleus. We also compared the magnitudes with the SDSS 
DR5 magnitudes. All quasars, except F2M0729-I-3336 
had SDSS imaging. Whenever the HST and SDSS mag- 
nitudes differ greatly, it means that the SDSS wasn't able 
to separate different components from the actual quasar 
due to its poor resolution. Table [3] presents the magni- 
tudes. 

The astrometry of ACS often has errors of up to 1.5 
arcseconds, so we shifted the WCS of the images us- 
ing SDSS stars in the field as references. This greatly 
improved our astrometry, to within a radial RMS er- 
ror of 0.1 arcseconds. We then compared these newly 
calibrated images with the FIRST data to confirm that 
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TABLE 3 

Comparison of total magnitudes 







mag 


'la 




mag / 










mag 








magi 






magg 






Source 


Aperture 


(3") 


Aperture 


(3") 


Isophot 




Isophot 




SDSS (petro) 


SDSS (petro) 


F2M07294 


-3336 


19.17 


± 


0.05 


21.71 




0.15 


19.20 ± 





05 


21.82 


± 





15 














F2M08254 


-4716 


21.04 




0.11 


22.45 




0.21 


21.05 ± 





11 


22.50 


± 





21 


18.93 ± 





08 


21.12 ± 





17 


F2M08304 


-3759 


18.45 




0.03 


20.24 


± 


0.07 


18.38 ± 





03 


20.19 


± 





07 


18.33 ± 





03 


20.06 ± 





04 


F2M08344 


-3506 


18.57 


± 


0.04 


20.31 


± 


0.08 


18.57 ± 





04 


20.39 


± 





08 


19.07 ± 





03 


20.98 ± 





08 


F2M08414 


-3604" 


20.45 


± 


0.09 


24.06 


± 


0.44 


20.23 ± 





08 


25.53 


± 





51 


20.10 ± 





09 


21.98 ± 





23 


F2M0915+2418 


19.53 


± 


0.06 


20.68 


± 


0.09 


19.56 ± 





06 


20.77 


± 





09 


19.95 ± 





06 


20.63 ± 





04 


F2M10124 


-2825 


20.23 


± 


0.08 


22.67 


± 


0.23 


20.24 ± 





08 


22.77 


± 





24 


20.38 ± 





63 


22.66 ± 





64 


F2M11134 


-1244 


18.60 


± 


0.04 


21.02 


± 


0.11 


18.64 ± 





04 


21.03 


± 





11 


18.78 ± 





03 


20.78 ± 





07 


F2M1118- 


-0033 


19.46 


± 


0.06 


21.71 




0.15 


19.31 ± 





05 


21.69 


± 





15 


19.41 ± 





06 


21.49 ± 





13 


F2M11514 


-5359 


20.05 


± 


0.07 


21.41 


± 


0.13 


20.05 ± 





07 


21.50 


± 





13 


20.06 ± 





11 


21.31 ± 





10 


F2M15074 


-3129 


19.48 


± 


0.06 


21.73 


± 


0.15 


19.39 ± 





05 


21.80 


± 





15 


19.73 ± 





12 


21.66 ± 





13 


F2M15324 


-2415'' 


19.46 


± 


0.06 


21.94 




0.16 


19.20 ± 





05 


22.06 


± 





17 


19.09 ± 





06 


20.78 ± 





10 


F2M16564 


-3821 


20.72 


± 


0.10 


23.65 


± 


0.36 


20.69 ± 





10 


23.68 


± 





36 


20.76 ± 





20 


23.54 ± 





92 



Note. — ^ HST magnitude quoted for brightest component. The different components are not resolved in SDSS. 
HST g -band magnitute quoted for brightest component. The different components are not resolved in SDSS or HST 
I-band. 



TABLE 4 
Sources with multiple components 



Source Comp. Separation 



Comments 



0841+3604 

1012+2825 
1118-0033 
1532+2415 
1656+3821 



1.70"/9 

0.15"/1 
0.40" /2 
0.75"/4 
0.80" /5 



2 kpc 
,2 kpc 

8 kpc 

9 kpc 
,8 kpc 



Wide separation 
Resolved in g -band 
One component stellar 
One component stellar 
2 components stellar 



the radio emission comes from the central nucleus or, 
in some cases, locate the likely position of the nucleus. 
The FIRST and SDSS surveys are fit to the same as- 
trometric system and both have position errors less 
than 0.2 arcseconds. In most cases the radio emission 
is associated with the brightest central optical compo- 
nent. For F2M1656-I-3821, the radio emission comes 
from the brightest of the three components. Yet for 
F2M0841+3604 the radio emission emanates from the 
center between the two brightest points. The lower a 
contours are skewed a little towards the brightest source. 
It is not clear where the radio emission is coming from 
in this source. 

4. ANALYSIS AND RESULTS 

4.1. Properties of the Quasars 

We performed quasar point-source subtraction and 
host galaxy fitting for the red quasars. We did per- 
form the subtraction for the most irregular merging sys- 
tems in our sample (F2M0841+3604, F2M1113-H1244, 
F2M1532+2415, F2M1656-f 3821); however, these 
sources have multiple components or are extremely ir- 
regular, so the host galaxy fitting will likely be wrong. 
For the fitting and subtracting we used our own IDL pro- 
gram fithost. I t is b ased on the 2D fitting technique of 
iMcLure et al.l (|2000( l and has been successfu lly applied 
to ground-based AO data (jLacv et al.l l2002\ The pro- 
gram simultaneously fits and subtracts both PSF and 
host galaxy; the host galaxy fitting is described in sec- 



tion [421 

PSF stars were observed at the end of the HST orbits, 
immediately after the quasars. Stars were chosen to be 
close to the quasars, so as to avoid extra overhead due to 
guide star acquisition, and to be bright enough to have 
a similar fluence to the quasar in short exposures. Ob- 
servations were made with the stars as close as possible 
to the same position as the quasars on the detector, and 
with the same dither offsets. The only drawback to this 
technique was that the ACS IK subarray had to be used 
for both the quasars and the guide stars to avoid exces- 
sive readout overhead. Since we were mostly interested 
in the quasars and their close companions this was not a 
major problem. 

The fithost program was used in a mode where it au- 
tomatically tries to fit the PSF-function to the brightest 
point in the image and then scale and remove it, so that 
the host galaxy can be subsequently fit. The center of 
the quasar and the host galaxy need not necessarily co- 
incide. Starting parameters for the fithost program were 
derived using the IRAF program imexam and then re- 
fined with a second run. The first 2 postage stamp im- 
ages in Figure [3] show the red quasars before and after 
PSF-subtraction. For the sake of brevity we chose to 
only show the /c-band results, since they have most con- 
tribution from the quasar and most incident flux. In the 
Ic-band the flux from the PSF is typically similar to that 
of the host galaxies, but concentrated in the center of the 
image, which is why in Figure [3] we don't see glaring dif- 
ferences in the PSF-subtr acted images from the original. 
In the g -band, the quasar is even fainter. The magni- 
tude of the PSF/quasar after fitting is found in the first 
two columns of Table [5] 

The first result we can derive from the PSF fitting of 
the HST images is that the quasars are actually more 
obscured than deduced from the spectrum. In the HST 
g -band, which corresponds to rest-frame UV, we almost 
only see host galaxy contribution with the quasar being 
almost entirely extinguished. The g -Ic colors the fitted 
quasar PSFs are also redder than for the total system. 
While the mean SDSS g'-i' color for the total system is 
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Fig. 3. — Jc-band results of PSF and Host Galaxy fitting. The first column shows an original postage stamp of the red quasar image, the 
second column the PSF-subtractcd image, the third the best-fit elliptical model, the fourth the residual after subtra<;ting the model and 
the fifth the radial surface brightness profile with the solid lines representing the elliptical fit. 
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Fig. 3. — cont. 
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Fig. 3. — cont. 
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1.92, the average g -Ic color for the quasar (just the PSF) 
is 0.61 magnitudes redder at 2.53. 

We then calculated the reddening for only the quasar 
by reddening the slope of the FBQS composite (account- 
ing for emission lines that might have fallen within the 
passbands) to fit the g -Ic colors made by the PSF us- 
ing the reddening curve described in section 12.11 Most 
of the objects which already had large E{B — V) red- 
denings around 1.0 increase their reddenings and colors, 
while the objects where the reddening of the total system 
were close to 0.5 on average don't increase or decrease 
their reddening. Interestingly, the objects for which the 
fit of the reddened quasar composite broke down, mostly 
decrease their reddenings, implying that the host galaxy 
itself has enough starlight that is responsible for some of 
the red colors of our objects. There could also be scat- 
tering of the nuclear light into the host galaxy itself, but 
we assume that the deduced PSF magnitude is the total 
quasar magnitude, since we can't assess the amount of 
scattering from the optical data alone. 

Figure m shows the shifts of the objects colors and red- 
denings. The total magnitudes and reddenings are rep- 
resented by the filled circles and the stars represent the 
quasar/PSF. Two of the quasars (F2M0841+3604 and 
F2M1656-f3821) were not resolved in the SDSS, so the 



PSF magnitudes are much fainter than implied by SDSS. 
We mark this objects in Figure H] in green color, their 
shifts are likely to be wrong, but we include them in the 
Figure for completeness. 

From the g -Ic colors, magnitudes and reddenings, we 
can then calculate the true luminosity for the quasar, 
since the absolute magnitudes calculated in section 12.11 
from the spectrum had contamination from the host 
galaxy, especially in the bluer passband. Table [5] gives 
the new luminosities (corrected for reddening) deter- 
mined from the quasar colors and magnitudes. The 
quasars are a bit more luminous on average, due to the 
fact that the quasar is more obscured than had been 
deduced from the spectrum. The quasars on average 
are still well above the quasar /Seyfert divide and much 
more luminous than the red quasar study conducted by 
iMarble et al] (p003l ). so our claim that we are missing a 
large population of obscured AGN is still valid. 



4.2. Properties of the Host Galaxies 

Before modelling the full host/quasar system, we ob- 
tained a constraint on the host galaxy magnitude. We 
did this by constraining the scaled PSF subtraction so 
that the residual was approximately flat in the center 
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TABLE 5 
Quasar (PSF) properties 



Source 




PSF map 


nitudcs 








g-I colors 


E{B - 


V) 






Luminosities (Mb) 






Ha 






Maf 








Total 


PSF 


Total 


PSF 






Total 


PSF 






















(SDSS) 


(HST) 


(spectrum) 


(HST) 




spectrum 


HST 


F2M0729+3336 


20 


15 


± 





07 


25.20 







73 


2.62" 


5.05 


0.83 ± 0.22 


1.18 ± 





16 


-24.96 


-26.83 


F2M0825+4716 


22 


33 


1 

it 


U 


zi 


23.65 







36 


O 1 o 

Z. iy 


i.OZ 


0.50 ± 0.35 


0.40 ± 





14 




01 no 




19 


51 


± 





06 


22 62 


-L 


n 
u 


22 


1.73 


3.11 


n sn -1- n 1 ^ 
u.ou in u.iu 


1.29 ± 





12 


-22.73 


-23.07 


F2M0834+3506 


18 


95 


± 





04 


21.17 


± 





11 


1.91 


2.22 


0.53 ± 0.10 


0.87 ± 





06 


-22.90 


-23.80 


F2M0841+3604 


21 


88 


± 





17 


25.11 


± 





70 


1.88 


3.23 


1.05 ± 0.64 


1.10 ± 





27 


-21.03 


-21.61'' 


F2M0915+2418 


19 


88 


± 





07 


20.99 


± 





10 


0.68 


1.11 


0.53 ± 0.36 


0.31 ± 





04 


-23.69 


-24.43 


F2M1012+2825 


21 


17 


± 





12 


23.50 


± 





33 


2.28 


2.33 


0.66 ± 0.12 


0.60 ± 





12 


-24.37 


-24.26'' 


F2M1113+1244 


19 


11 


± 





05 


22.87 


± 





25 


2.00 


3.76 


1.01 ± 0.24 


1.17 ± 





09 


-25.08 


-25.45 


F2M1118-0033 


22 


DO 


± 





18 


23.33 


± 





32 


2.08 


1.33 


0.72 ± 0.26 


0.48 ± 





15 


-22.97 


-21.76 


F2M1151+5359 


20 


48 


± 





09 


22.12 


± 





18 


1.25 


1.59 


0.47 ± 0.13 


0.47 ± 





07 


-23.31 


-23.78 


F2M1507+3129 


19 


78 


± 





07 


22.67 


± 





23 


1.93 


2.89 


0.54 ± 0.13 


0.69 ± 





08 


-24.27 


-26.24 


F2M1532+2415 


21 


14 


± 





12 


24.42 







51 


1.69 


3.28 


0.90 ± 0.53 


1.12 ± 





21 


-20.46 


-22.45'' 


F2M1656+3821 


23 


46 


± 





35 


25.11 


± 





70 


2.78 


1.65 


0.55 ± 0.57 


0.54 ± 





28 


-20.86 


-20.66'' 



Note. — "^No SDSS photometry, wc use HST values. ^ Luminosity is only for one nueleus. 



The position of the PSF and the galaxy nucleus was not 
held fixed; so in total we fit the flux, angle, axial ratio 
and position of the host galaxy to the image. We tried 
to fit both elliptical and exponential profiles to the im- 
ages, however in all, but one case (F2M1 118— 0033), the 
elliptical profile was best, so all our quoted values are for 
an elliptical profile. 

Figure[3]shows the fitted elliptical galaxy and the resid- 
uals after PSF and model-galaxy subtraction in the third 
and fourth panels respectively. Note that none of the 
galaxies are a perfect fit to an elliptical, however, there 
were some which came quite close, like F2M0915-I-2418. 
Also, two of the host galaxies, while they did not fit 
the elliptical profile well, did not show signs for inter- 
action (F2M0834-H3506 and F2M1151-h5359). The frac- 
tion of 11 out of 13 systems showing interaction is much 
higher than the 2M ASS red AGN fraction of one third 
([Marble et al.l[2003h and other unobscured quasar host 
galaxy studies. The fifth panel in Figure [3] shows the 
radial surface brightness profiles for the PSF subtracted 
images and the solid line is the model applied with which 
they were fitted. 

In the following paragraphs, when we discuss the host 
galaxies we mean the PSF-subtracted images and not 
the modeled host galaxies, to assess their magnitudes, 
morphologies and comparison the the quasar. 

We measured the magnitudes of the host galaxies 
with Sextractor following the conventions above. The 
isophotal-magnitudes are quoted in Table [6l They were 
calculated in the same way the magnitudes were calcu- 
lated in Section [3l Note that in the cases where the total 
magnitude of the host galaxy is fainter than the limit de- 
noted by the "zero" magnitude, the total host magnitude 
only refers to the central component of the galaxy, as the 
host galaxies often showed so large interactions that they 
are broken up and have various components. 

The colors of the host galaxies are not particularly 
blue, so the host galaxies do not show large amounts 
of unobscured star formation. However, some of the ex- 
tended components are blue, indicating the presence of 
either merger-induced star formation, or light associated 



Average shift: 

B -I,: 0.61 
E(B-V): 0.085 



ii 



/? 



0.6 0.8 
E(B-V) 



Fig. 4. — g -Ic color vs. E(B — V) reddening. The filled dots 
represent the magnitudes and reddenings for the total system. The 
magnitudes are from SDSS data and don't differ much from HST 
photometry and the E(B — V) reddening is from spectral fitting. 
The open stars represent quasar (PSF) color and reddenings de- 
rived from HST PSF photometry. We included lines to show how 
the values shifted. On average the quasar has a 0.61 mag redder 
color and higher E{B — V) of 0.085 mag. 



of the quasar host within a 3 pixel radius and declined 
monotonically outside that radius. We also obtained the 
lower limit of the host galaxy magnitude by subtracting 
the PSF until the central pixel was zero. The "mono" 
and "zero" magnitude are found in Table [G] 

The host galaxies were modeled by fitting PSF plus 
galaxy model profiles (convolved by the PSF) by mini- 
mizing . Close to the center of the source, systematic 
errors from the PSF subtraction dominate, so to prevent 
those errors to dominate the fit, the inner 3 pixel radius 
was downweighted by a factor of 0.5 to ensure that the 
surface was fairly uniform across the fitting aperture. 
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E2M{B 15+2418 K/II = 1 



F2M09 15+24 IS K/H = 3 



F2M(,)9 15+24 IS N/H = K.t 




R2MI11SK)033 N/H=l 
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F2M 1 532+24 15 N/H = I F2M 1 532+24 15 N/H = 3 



F2M 1 532+24 1 8 N/H = 10 






Fig. 5. — Host galaxies with artificial PSF added. The images show different nucleus to host ratios in increasing order. When the nucleus 
has the same magnitude as the host galaxy, which is the case for the red quasars in this sample, the interactions are clearly visible, such 
as the tidal connection in F2M1118— 0033 or the strong halo in F2M1532+2415. With a nucleus to host ratio of 10, most of the features 
have disappeared and the PSF clearly dominates, as is typical in unobscured quasar sampled. 



with the quasar. This quasar light could be either scat- 
tered continuum, or extended Une emission. The lack of 
a clear ionization cone morphology (with the possible ex- 
ception of F2M08304-3759), however, would suggest that 
most of this blue light is from young stars. 

We also measured the nucleus to host (N/H) ratios 
from these magnitudes (TablelHlfor the /c-band N/H). As 
discussed in section ITTl the N/H for the g -band will be 
lower, because the quasar is more extinguished and most 
of the light in that band will come from the host. Overall, 
the luminosity of the hosts are on average, a little larger 
than those of the nuclei, such as the sample of IR-excess 
quasar sample of Suracc, Sanders & Evans (2001). The 
relatively low N/H are in sharp contrast to the high N/H 
ratios of the luminous quasar sample of iMcLure et al.l 



(|1999l ). which is tied to the iDunlop et all (|2003D sam- 
ple. However the N/H ratios of our quas ars are gener- 
ally h igher than the red quasar sample of iMarble et al.l 
(|200l. 

We then calculated the absolute magnitudes {Mb) of 
the hosts adopting a 1 G yr post starburst model from 
iBruzual fc CharlotI (|2003D for the K-correction. At the 
typical redshifts o f our sample (z ~ 0.7) Mg = —20.30 
(jllbert et al.ll200^ . so the host galaxies are around or a 
bit above L* luminosity. 

One measure of t he extent of the host galaxy is the 
Petrosian radius Vp (jPetrosianl Il976f ) . Very irregular or 
interacting systems tend to have high Petrosian radii, 
while compact ellipticals have a low r^. The Petrosian 
radius is defined when the ratio of he surface brightness 
at the Petrosian radius and the average surface bright- 
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ness at radii below reach a certain value ij: 



Kr < rp) 



Following SDSS conventions, we set 7] = 0.2. Ta- 
ble [6] quotes the Petrosian radii of our systems. With 
one exception (the very compact F2M0729-I-3336), most 
have Petrosian radii around l"-2.5", which are normal 
for galaxies at those redshifts. When comparing those 
numbers to the surface density plots in Figure [31 one 
can notice that the nucleus (r < r^) of the host galaxy 
often has an elliptical profile; only the most irregular sys- 
tems such as F2M0841-I-3604 are irregular at low radii. 
Mostly, the interaction features only appear well beyond 
the Petrosian radius and usually have low surface bright- 
ness. 

This could be an indication why many authors con- 
clude that the host galaxies of luminous quasars are 
fit be st mostly by elliptical profiles (e.g. (jFlovd et al.l 
|2004[ )). but find so few merger remnants. Perhaps with 
the overpowering AGN (high N/H ratios), the low sur- 
face brightness features or interaction morphology signa- 
tures are lost. Recent ly, very deep imaging of some of 
the quasars studied bv lDunlop et al.l (j2003 ). show either 
low surface brightness tida l tails or other merg er rem- 
nants in the form of shells (|Bennert et al.|[2b0 6') further 
supporting this point. The red quasars in this study 
are quite extinguished, so the host galaxy features are 
easier to discern. We tested how the tidal features can 
disappear in the presence of a bright quasar by adding 
point sources of different brightness (1,3,10 N/H ratios) 
to the host galaxy image. When a PSF with a N/H ratio 
of 10 is added, the quasar is so bright that in all cases 
with exception of F2M0841-I-3604 the red quasars lose 
most of their interaction features and the quasar dom- 
inates the image (Figure O. Even so, after performing 
PSF subtraction on a N/H=10 quasar -I- host system, the 
interactions should reappear. Also the red quasar host 
galaxies are in an earlier stage of the merger as we will 
see next. 

As we already commented in section |3] the fraction of 
our red quasars showing interaction in their host galax- 
ies is very high. Furthermore, all of the dust reddened 
quasars, that is the quasars where the dust reddening 
template fit best, show interaction. E(B-V) seems to cor- 
relate weakly with the amount of interaction; the more 
obscured the quasar is, the more disturbed the morphol- 
ogy of the host galaxy is. We can arrive at these conclu- 
sions more or less "by eye" , but we support our claims 
by parameterizing the host galaxy morphology with their 
Gini coefficients and their Concentration indices. 



The Gini coefficient is a non-parametric approach to 
classifying a galaxy, and can therefore be easily ap- 
plied to irregular galaxies or highly merging systems. 
It is a measure of the cumulative distribution of a 
galaxy's pixel values and is a good alternative approach 
to quantifying the amount of interactions in galaxies. 
The Gini coefficient is correlated with the concentra- 
tion index of galaxies for spiral and elliptical galax- 
ies ({Abraham, van den Bergh fc Naiil [200 3i) and anti- 
correlated with the Concentration Index for ULIRGs and 



other highly irregular systems (jLotz. Primack &: Madaul 
l2004f ). A high Gini coefficient therefore indicates either 
a high level of interactions or a highly concentrated ellip- 
tical galaxy, while a low Gini coefficient indicates either 
spiral galaxies or low surface brightness galaxies. 

We calculated the Gini coefficient of the host galax- 
ies of the red quasars by follow i ng the conventions in 
[Abraham, van den Bergh &: Naiil (|2003( ). Based on the 
formalism bv lGlasser et afT i 19621 ) . if we sort the pixels 
Xi's flux into increasing order, the Gini coefficient can 
be calculated by: 



G 



Xn{n 



1 " 
1-1) ^ 



1)X, 



(3) 



However, we took notice of the warning by 



[Lotz. Primack &: Madaul ([2004D and tried to create seg- 
mentation maps at the fi{rp), the flux threshold above 
which pixels are assigned to the galaxy. For that we 
cut out galaxy postage stamps cutouts by eye and in- 
cluded only pixels within that cutout that had a surface 
brightness higher than n{rp). The Gini coefficient was 
then computed of the distribution of absolute flux values 
of the pixels which corrects for the areas in the center 
which were PSF-oversubtracted. The Gini coefficients G 
can be found in Table [6l 

The Concentration Index measures the ratio of the flux 
withing an inner radius to that within and outer circular 
or elliptical aperture. The main difference in the defi- 
nitions by different authors is the choic e of the radii or 
semimajor axes of the two apertures. [Conselice et al.l 
(|2003f l adopts a ratio between th e radii containing cer - 
tain percentages of the total light. [Abraham et al.l ( 19941 ) 
uses a flux ratio between two normalized radii a and 1, 
where E(l) is the area encompassing 2a flux and a is 
typically set to 0.3. 



C = 



EE 



i,jeE{a) ^ij 



E Ei,j 



(4) 



Since [Abraham, van den Bergh fc NaiH (|200l has al- 
ready linked the Gini coefficient to the Concentration 
index via a unity slope, we decided to use the Con- 
centration inde x defin ed in equation [31 instead of the 
[Conselice et al.[ ([20031) definition. We chose the total 
flux as the flux within l.Sr^ and E(a) as the flux within 
0.45rp. From that ratio we are able to obtain C, which 
is quoted in Table [SI 

Figure [SI shows the Gini coefficient plotted against 
the E{B — V) reddening derived from the HST quasar 
magnitudes on the left side and the Concentration In- 
dex over the same E{B — V) on the right side. While 
there is a weak correlation between the Gini coefficient 
and the reddening, we find no correlation between the 
Concentration Index and E{B — V). As mentioned be- 
fore, high Gini coefficient and low Concentration indices 
are an indication for interaction or merger, since the 
brightest flux pixels are highly concentrated in a few 
pixels (high G), but those pixels are not in the cen- 
ter (low C). This seems to be the case in the systems 
with the highest E{B — V) in our sample, while the 
hosts of quasars with lower E{B — V) tend to have 
both Gini coefficients and Concentration indices more 
consistent with those of normal, undisturbed galaxies 
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TABLE 6 
Host galaxy properties 



oourcG 




Host magj^ 




Host mag / 


Luminosities 




Morphologies 






zero 


mono 


total 


zero 


mono 


total 


N/H" 


Mb 


L* 


fnix) 


G 


C 




20.71 


19.97 


19.78 ± 0.06 


21.80 


21.80 


21.78 ± 0.15 


0.711 


-22.32 


6.4 


03.8 


0.69 


0.36 


F2M0825+4716 


20.07 


19.99 


21.48 ± 0.14 


21.71 


21.58 


22.95 ± 0.26 


0.470 


-20.02 


0.8 


29.3 


0.54 


0.42 


F2M0830+3759 


18.96 


18.75 


18.86 ± 0.04 


20.53 


20.22 


20.31 ± 0.08 


0.550 


-21.02 


1.9 


25.9 


0.57 


0.41 


F2M0834+3506 


19.96 


19.87 


19.77 ± 0.06 


22.66 


22.56 


21.13 ± 0.11 


2.089 


-20.42 


1.1 


37.5 


0.52 


0.55 


F2M0841+3604 


19.96 


19.90 


20.40 ± 0.09 


21.57 


21.54 


23.04 ± 0.27 


0.256 


-20.22 


0.9 


53.7 


0.52 


0.30 


F2M0915+2418 


20.72 


20.64 


20.82 ± 0.10 


25.58 


23.83 


22.91 ± 0.25 


2.355 


-20.90 


1.7 


24.0 


0.49 


0.38 


F2M1012+2825 


20.56 


20.44 


20.87 ± 0.11 


22.12 


21.99 


23.49 ± 0.33 


0.759 


-21.15 


2.2 


19.8 


0.59 


0.52 


F2M1113+1244 


20.06 


19.76 


19.75 ± 0.06 


21.12 


21.87 


21.31 ± 0.12 


1.803 


-21.28 


2.5 


18.6 


0.68 


0.55 


F2M1118-0033 


19.34 


19.24 


19.51 ± 0.06 


20.95 


20.91 


22.70 ± 0.22 


0.101 


-21.53 


3.1 


22.6 


0.59 


0.42 


F2M1151+5359 


20.03 


19.99 


20.99 ± 0.11 


22.70 


22.46 


22.43 ± 0.20 


1.528 


-20.41 


1.1 


49.6 


0.39 


0.41 


F2M1507+3129 


20.10 


20.07 


20.96 ± 0.11 


22.24 


22.09 


22.45 ± 0.21 


2.858 


-21.18 


2.2 


48.5 


0.57 


0.49 


F2M1532+2415 


19.15 


19.05 


19.39 ± 0.05 


20.77 


20.69 


22.20 ± 0.18 


0.194 


-21.27 


2.4 


33.4 


0.60 


0.40 


F2M1656+3821 


20.58 


20.45 


20.77 ± 0.10 


22.22 


22.14 


24.96 ± 0.62 


0.039 


-20.47 


1.1 


20.3 


0.59 


0.44 



Note. — ^Nucleus to host ratios from /c-Band magnitudes. 
Zero and monoton are model magnitudes, therefore they don't have errors. 



lAbraham. van den Bergh fc Naiil (|2003[ ). On the lower 
panel of Figure [6] the Concentration Index is plotted 
against the Gini coefRcient, with the unity slope shown. 
Only the two undisturbed galaxies (F2M0834+3506 and 
F2M115H-5359) are above this slope, with the rest hav- 
ing high Gini coefficients in relation to their Concentra- 
tion Indices, again implying high interaction features. 

While we could already see those results "by eye", 
these galaxy parameters are affirmation to the result that 
the higher the reddening in the quasar, the higher the 
chance of interaction in the quasar host. Furthermore, 
if the reddening in the spectrum fits a dust reddening 
template well, it is also a good indicator for evidence of 
interaction. Nonetheless, these results are based only on 
13 highly heterogenous systems, so further observations 
are needed to improve number statistics and to confirm 
this claim. 

5. NOTES ON INDIVIDUAL QUASARS 

While most of the quasar host galaxies show some de- 
gree of irregularity and evidence for merger, it is worth 
it to look at the quasars on an individual basis to further 
inspect some peculiarities in some of the objects. 

5.1. F2M0729+3336 

This is the only quasar for which SDSS information was 
not available. Even though we detected only 6 photons in 
the X-ray observations, all of them are in the hard band, 
giving this source a hardness ratio of 1.0 and therefore 
implying a high column density. 

The spectrum has strong Ca H-fK absorption lines, 
but the spectral reddening does fit a dust extinction law 
quite well longward of about 5000 A , so the host galaxy 
contribution to the red optical light in the spectrum is 
not going to make a large impact. 

This particular system has a very compact elliptical 
host, with a very small Petrosian radius and large Gini 
coefficient. Beyond that however, large tidal tails extend 
out to over 20 kpc on both sides of the host, which is 
evidence for a recent merger event. This makes the to- 
tal size of the host galaxy over 35 kpc. In the g -band, 
the quasar (PSF) contribution is almost non-detectable. 



making this the system with the largest shift in g -Ic 
color of the PSF relative to the total emission from the 
host. 

5.2. F2M0825+4716 

This looks like a Type-2 quasar in the optical spec- 
trum. The color of the quasar (PSF) would suggest that 
also. However, the Spex IR spectrum shows broad quasar 
lines, however and a much redder spectral slope than im- 
plied by the optical continuum, which is probably dom- 
inated by host galaxy emission. The image shows that 
while the central component is fit well by an elliptical 
galaxy, there is a tidal bridge linking the central galaxy 
to a companion nearby. 

The optical spectrum shows double-peaked narrow 
emission lines, separated by about 600km s^^. The 
higher redshift set has relatively strong Balmer lines and 
[Oii]3727 emission at z = 0.803, and probably corre- 
sponds to the systemic redshift of the host galaxy. The 
blueshifted component at z = 0.800 is of higher ioniza- 
tion and probably corresponds to a high ionization out- 
flow. 

5.3. F2M0830+3759 

This quasar has the lowest redshift in our sample, 
therefore features close to the nucleus will be more easily 
observable. The Hubble image in fact shows a lot of ir- 
regularity near the nucleus with shells of material along 
the major axis, either side of the nucleus. They almost 
cancel each other out in the radial profile plot, resulting 
in a deceptively smooth radial profile, while in reality the 
host galaxy is very disturbed. 

F2M0830-I-3759 was also obse rved with Chand ra and 
has a very high X-ray flux (Ur rutia et al.l [20051 ). The 
X-ray spectrum shows a very broad Iron Ka line, hint- 
ing that we are looking into the broad line region of 
the quasar. It has a moderate absorption, but a higher 
than usual spectral slope (F « 2.9). This is also the 
only object for which we have a gas:dust ratio, which 
is about 5 times Galactic value. While dusty, it is well 
within the 3-100 time s Galactic value range deduced by 
iMaiolino et"all (|200lD . 
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Fig. 6. — Upper left graphic: Gini coefficient vs. reddening. There is a weak correlation between the Gini coefficient and the reddening 
(dashed line). High Gini coefficient indicate either highly irregular systems, such as ULIRGs or very compact ellipticals. Upper right 
graphic: Concentration index vs. reddening. For single systems the Gini coefficient and the Concentration index are correlated. The fact 
that for high reddening we don't necessarily have high Concentration indices means that the systems are highly disturbed, usually with 
more than one nucleus. Lower graphic: Concentration index vs. Gini coefficient. No correlation is found in our sample in contrast to other 
studios of single galaxies, indicating that a many of the red quasar hosts are in a merger status. The dotted line is the unity correlation 
found by .Abraham, van den Bergh Nair (,2003i ) . All but two host galaxies (the two undisturbed ones) of the F2M red quasars have 
higher Gini coefficients than implied from their Concentration Indices. 



5.4. F2M08S4+3506 

The red quasar system has ahiiost no evidence for in- 
teraction. The profile could fit an addition of elliptical 
profiles, but the morphology could also represent merger 
remnant shells. However, there are no extended features 
or tidal tails. 

It is unclear whether the blue component just to the 
North of the quasar is associated with it, or whether it is 
an unrelated object close to the line of sight. There are 
broad Ha and H/3 emission lines ^^5000 km s~^ blueward 
of the quasar's redshift, which could be associated with 
that component and/or with the quasar. 

5.5. F2M0841+3604 

This object shows the most irregular morphology in 
the host galaxy of the Hubble images. While the optical 



image clearly shows two bright nuclei, of which either or 
both could be the quasar, the radio images (from FIRST 
and also 6cm VLA A-array observations) show a steep 
spectrum point source right in between the two bright 
optical components. This raises the question of where 
the active nucleus actually is. 

In addition to the HST and VLA images, 
F2M0841-h36()4 ha s been observed with Chandra 
(jUrrutia et al.l |2005| ) . The Chandra photons are very 
hard. Unfortunately, the X-ray image had only 7 counts, 
so the position error on where the X-ray emission comes 
from is quite large. However, the Chandra photons ap- 
pear to be spread out, so the X-ray emission could come 
from both nuclei. Overall, this is our most spectacular 
example of a young evolutionary merger state. 

5.6. F2M0915+2418 
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F2M0915+2418 is the quasar which fits an elliptical 
profile best. However, a clear tail is seen extending to th 
East indication a recent interaction. 

5.7. F2M1012+2825 

This red quasar displays two nuclei only 0.15" (1.2 kpc) 
apart, which would have been missed in ground-based 
observations. The two nuclei are point-like and espe- 
cially easy to detect in the g -band, but the nucleus ap- 
pears very "smeared" in the /c-band image. The fitting 
of quasar and host galaxy is therefore impossible. The 
quoted luminosities for the quasar are only for one com- 
ponent, so the total nuclear luminosity is much higher 
than Mb = -24.26. 

The rest of the host galaxy does not appear very dis- 
turbed and also here the central component fits an ellip- 
tical profile, were it not for that second nucleus. There 
is some blue emission north of the galaxy, but again, we 
do not know if it is associated with the quasar host. 

5.8. F2M1113+1244 

While the position of the quasar is quite certain, be- 
cause of the clear stellar-like feature, it is noteworthy 
that the rest of the host galaxy consists of several bright 
knots, with a low surface brightness tidal tail extending 
to 30-35 kpc from the galaxy. 

5.9. F2M1118-0033 

The image for this red quasar shows it to be a clear 
merger with a nucleus and two compact components, one 
near the nucleus and the other one at the end of the two 
tidal tails trailing the system. This is the only system 
for which an exponential profile fit the host galaxy better 
than an elliptical (but only near the nucleus). 

The optical spectrum for F2M1118-0033 shows a 4000 
A Balmer break, and a red spectral slope which appears 
to be dominated by the hos t galaxy. The infrared spec- 
trum iGlikman et al.l (|2004l ) also fails to show obvious 
broad emission lines, but does show a significant rise in 
K-band, probably from quasar light. 

5.10. F2M1151+5359 

F2M1151-f 5359 is the quasar also fits an elliptical pro- 
file quite well, but there is still some residual in the host 
galaxy after model subtraction in the /c-band. The blue 
band fits almost perfectly. So this would be comparable 
to other "normal", bright, blue quasars, which have an 
undisturbed elliptical galaxy as a host. The spectrum 
shows asymmetric profiles in the emission lines and nar- 
row Mgll absorption lines very close to the quasar red- 
shift. 

5.11. F2M1507+3129 

This is the highest redshift red quasar imaged with 
ACS. The spectrum fits dust reddening quite well. Also 
here the central component can be well fit by an el- 
liptical galaxy, but there is a low surface brightness 
tidal bridge connecting it to the red companion nearby. 
F2M1507-f 3129 is thus a very large (35 kpc) and spec- 
tacular merger at z~l. 

5.12. F2M1532+2415 



HST imaging of F2M1532-I-2415 displays a truly dis- 
turbed host, with many components and no clear ev- 
idence of where the quasar is located among the con- 
stituents of the host galaxy. Radio maps show the syn- 
chrotron radiation coming from the central red com- 
ponent. They also show that this quasar is actually 
a classical Fanaroff-Ril ey Type II radio double (FRII, 
iFanaroff fc RilevI ()1974l )'). with the outer components be- 
ing over one arminute away from the central source. 

The optical spectrum shows a very narrow II/3 line and 
no broad Mgll line, but in the infrared Pa/3 is broad. The 
quasar host shows extreme interaction. 

5.13. F2M1656+3821 

This optically-faint object shows almost only host 
galaxy emission in the optical part of the spectrum, but 
in the near-infrared there is a very red continuum and a 
broad Pa/3 line. There are three components in the HST 
image, two point-like and one elliptical. The radio flux 
points to the emission coming from the brightest (SW) 
optical point source. The morphology of F2M1656-|~3821 
could suggest that the system is a gravitational lens with 
the lensing galaxy lying between the two point sources. 
The reddening for this system might stem from the lens- 
ing galaxy and not the quasar itself, as is implied by 
the host galaxy's very red colors. However, the lack of 
any arc-like morphology in the outer optical components 
and the lack of spectral features at lower redshift in this 
object's spectrum argues against the lensing hypothesis. 

6. CONCLUSIONS 

We have observed and analyzed a sample of 13 
highly reddened, luminous quasars from the F2M survey 
'Glik man et all (|2007D with ACS in both Ic- and ^'-band 
to study their host galaxies. The images show interac- 
tions in 85% of the objects and clear evidence for mergers 
such as tidal tails and multiple nuclei. We fitted model 
PSF and galaxy profiles to the images. Within the nu- 
clear region, most galaxies fit an elliptical profile, in ac- 
cordan ce with quasar host studies such as lDunlop et al.l 
(|2003f ). but outside of the nuclear region this fit is no 
longer valid as merger features become clear. 

After performing PSF-fitting, the quasar displays even 
redder colors and reddenings than implied from SDSS 
imaging and ESI-spectra, showing that the optical mag- 
nitudes of the quasars are significantly contaminated 
by the host galaxies in low resolution data. Five of 
our galaxies show very dramatic, likely multiple merger 
morphologies, reminis cent o f the simulations of high- 
z quasar hosts by Li et al.l (|2006l ). There is a wide 
range in merger phase, with some of our quasar hosts 
in an apparently early part of the merger sequence 
where the merging galaxies still retain their identity 
(F2M0841-h3604), through to later stages with multi- 
ple nuclei (F2M1 113-1-1244; F2M1012-)-2825), and ob- 
jects at the final merger stage with only a single nucleus 
and a tidal tail remaining (F2M0915-f 2418). This would 
be consistent with quasar activity being triggered rela- 
tively early in an interaction, and continuing through the 
merger process. 

We calculated the Gini coefficients and the Concentra- 
tion indices of the host galaxies and found a correlation of 
the Gini coefficient with reddening. The systems with the 
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highest quasar E{B — V) did not have high Concentra- 
tion indices consistent with their high Gini coefficients, 
which imphes that the host galaxies have bright, compact 
components outside of the nuclear region. Using this as 
a measure of interaction, it seems that the amount of 
interaction is weakly correlated with how obscured the 
q uasar is. 

iFlovd et al.l (|2004[ ) studied a sample of 17 normal 
quasar hosts with HST which are closest to our sam- 
ple in terms of luminosity and redshift. Their sam- 
ple spanned —24 < My < —28 in quasar luminosity 
and 0.3 < z < 0.42 in redshift. Consistent with previ- 
ous studies, only a small fraction of interactions/mergers 
were seen (they have only one object out of 17 with def- 
inite signatures of an interaction). However, sensitive 
observations of the hosts of unreddened quasars of sim- 
ilar luminosity are difficult, so it is possible that some 
signs of interaction such as faint tidal tails and multiple 
nuclei close to the quasar will have been missed. Our 
correlation of reddening with morphology argues against 
this, but our observed correlation is only weak. However, 
Figure [5] also shows obvious signs of merger even at nu- 
cleus to host ratios of 10. Another possible source of bias 
is that, while none of the red quasars is radio-loud, these 
objects are all radio-selected, and fall into the "radio in- 
termediate" category. The fact that all of our quasars are 
disturbed in comparison to, e.g., the|Marble et al. (2003) 
sample, could be a selection effect. I White et al.i ()2006il 
find that quasars with redder colors than the SDSS com- 
posite tend to have higher radio-fluxes. 

These caveats aside, our results may explain why only 
about 30% of quasars with obviou s signs of mergers 
and interactions h ave been found (jCuvon et al.l 120061 : 
iMarble et al.l 120031 ) . Reddening by dust from the host 
galaxy seems to be the most likely explanation for the 
redness of the F2M quasars. The host morphologies and 
colors are consistent with dusty, merging galaxies. If 
these objects were red because the line of sight just grazes 
the torus we would expect less disturbed host morpholo- 
gies, in line with those of normal quasars. Host galaxy 
reddening of the quasar, whilst only mild in our cases 



(and much less than nuclear reddening by an edge-on 
torus) can nevertheless be sufficient to remove the objects 
from optically-selected quasar samples (or at best render 
them a small minority in such a sample). Given that 
these kind of dust-reddened Type-1 quasars like these 
make up about 30% of the total in mid-IR selected sam- 
ples, this means the fraction of quasar hosts associated 
with mergers has been significantly underestimated in 
the past. Our result is thus consistent with theories in 
which quasars start their lives obscured by dust, and only 
appear in optical or soft X-ray surveys after the dust 
along the line o f sight to the nucleus has been cleared 
by quasar winds (ISanders e t al."1989':'Sanders fc Mirabel 
flQOalSpriMel et al.ll2005nHopkins et al. 200(1) . 
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